Completeness of classical spin models and universal 
quantum computation 



Gemma De las Cuevas^, Wolfgang Diir^, Maarten Van den 
Nest§ and Hans J. Briegel^ 

Hnstitut fur Quantenoptik und Quanteninformation der Osterreichischen Akademie 
der Wisscnschaftcn, Innsbruck, Austria 

Hnstitut fur Thcorctische Physik, Universitat Innsbruck, Tcchnikcrstrafic 25, A-6020 
Innsbruck, Austria 

§ Max-Planck Institut fur Quantenoptik. Hans-Kopfcrmann-Str. 1, D-85748 Garching, 
Germany 

Abstract. We study mappings between distinct classical spin systems that leave the 
partition function invariant. As recently shown in [Phys. Rev. Lett. 100 110501 
(2008)], the partition function of the 2D square lattice Ising model in the presence of 
an inhomogencous magnetic field, can specialize to the partition function of any Ising 
system on an arbitrary graph. In this sense the 2D Ising model is said to be "complete" . 
However, in order to obtain the above result, the coupling strengths on the 2D lattice 
must assume complex values, and thus do not allow for a physical interpretation. 
Here we show how a complete model with real — and, hence, "physical" — couplings 
can be obtained if the 3D Ising model is considered. We furthermore show how to 
map general g-state systems with possibly many-body interactions to the 2D Ising 
model with complex parameters, and give completeness results for these models 
with real parameters. We also demonstrate that the computational overhead in 
these constructions is in all relevant cases polynomial. These results are proved by 
invoking a recently found cross-connection between statistical mechanics and quantum 
information theory, where partition functions are expressed as quantum mechanical 
amplitudes. Within this framework, there exists a natural correspondence between 
many-body quantum states that allow universal quantum computation via local 
measurements only, and complete classical spin systems. 
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1. Introduction 



Classical spin systems, such as the Ising and Potts models pQ, are widely studied 
in physics and appear in different contexts, ranging from magnetism to economic 
models [2]. Even in two spatial dimensions, these systems can show interesting and 
complex features — as is apparent in e.g. their phase diagrams — and they can be difficult 
to treat analytically or numerically [3]. 

One of the central characteristics of a spin system is its partition function, as many 
relevant macroscopic properties of the system can be derived from it. In this paper we 
concentrate on mappings between classical spin systems of possibly different size that 
leave the partition function invariant, and thus also the thermodynamical properties 
of the system. On the one hand, we will identify different models that have the same 
partition function and are hence in this sense equivalent. On the other hand, we will find 
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that there exist certain models that are complete. Colloquially speaking, a spin model 
Ai, defined on a certain class of lattices and exhibiting a certain range of coupling 
strengths, is said to be "complete" if the partition function of every other classical 
spin system Ai' can be re-expressed as a partition function of the model At, where 
the lattice size of M as well as its couplings depend on M> % Such completeness 
results were first established in [I], where it was shown that the 2D Ising model with 
inhomogeneous magnetic fields constitutes a complete model with respect to all q— state 
Potts-type models with pairwise interactions. That is, the partition function of every 
such Potts-type system (with arbitrary finite q, defined on an arbitrary graph and with 
arbitrary couplings) can be re-expressed as the partition function of an Ising model on 
a 2D square lattice which is polynomially larger than the original system, with suitably 
tuned parameters 0. 

The main tool to obtain these results (as first used in [I]) is a new connection 
between classical statistical mechanics and quantum information theory (see also [51 El 
H El El Cni [HI [121 H3l [HI HSl HEl HZl HE], where similar connections are studied). In ^ 
the partition function of a classical spin system is expressed as a quantum mechanical 
object, namely a quantum amplitude. More specifically, partition functions Z are 
shown to coincide with certain overlaps (inner products) between pairs of many-body 
quantum states: Z = (ot\ip). Here \ot) is a "product state", i.e. a state which does 
not contain any entanglement, whereas is a (typically highly entangled) "quantum 
stabilizer state" . Quantum overlaps of the above form appear naturally in a branch of 
quantum information theory called "measurement-based quantum computation" [19] 
(or "one-way quantum computation"). This a paradigm for quantum computation in 
which computations are carried out by performing sequences of (adaptive) single-qubit 
measurements on a highly entangled many-body "resource state". It is known that 
there exist certain stabilizer states, such as the two-dimensional (2D) cluster states [20] 
that are universal states for measurement-based quantum computation. That is, every 
quantum computation can be carried out by merely performing a suitable sequence of 
single-qubit measurements on such universal states. It is this universality feature of the 
2D cluster state, together with the aforementioned connection between stabilizer states 
and classical spin systems, that is used in [4] to establish the completeness of the 2D 
Ising model. 

In the present paper, this connection to quantum information theory will be further 
studied and exploited to strengthen the results of [I]. The completeness results presented 
in [1] exhibited some limitations, which will be here significantly relaxed. For example, 
in |4] the completeness of the 2D Ising model could only be obtained when allowing 
for negative- or complex-valued — instead of real and positive — Boltzmann weights, 
making it hence difficult to give a physical interpretation of this result. Furthermore, 

\ In this context we sometimes say that the partition function of the complete model "specializes" to 
the partition function of all other models. 

§ If the completeness only holds for a restricted set of models on certain graphs, it is specified by 
saying, e.g., "completeness for Ising models". 
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efficient reductions — i.e. requiring only a polynomial amount of operations in the system 
size — were in [I] only provided from Potts- type models to the 2D Ising model; for models 
which were not of Potts-type (for example, vertex models such as the six-vertex model, 
which involve four-body interactions), it was unclear whether such reductions to the 
2D Ising model could be made efficient as well. In the present paper we show that the 
above limitations can be overcome in several cases. The following is a summary of our 
main results. 

• The Ising model with inhomogeneous magnetic fields on a 2D square lattice (2D 
Ising model) is complete with complex parameters for all spin models. This includes 
(edge and vertex) models on arbitrary graphs (hence including lattices of arbitrary 
dimension) with q— level spins and k— body interactions, for every finite q and k. 
The overhead (i.e. the size of the 2D Ising model as compared to the original system) 
is polynomial in the number of spins, as long as q and k are bounded (i.e. arbitrary, 
but fixed). 

• The 3D Ising model without magnetic fields is complete with real parameters for 
all Ising models with magnetic fields on arbitrary graphs, and only a polynomial 
overhead is required. Similarly, we present a q— state edge model on a 3D square 
lattice which is complete with real parameters for all q'— state edge models with 
pairwise interactions as long as q' < q. 

• Going beyond Ising models as complete models, we extend the analysis to vertex 
models and we construct a certain 2D two-state vertex model which is complete 
for all models when allowing for complex parameters. When restricting to the real 
parameter regime, we construct a 2D q— state vertex model which is complete for 
all q'— state vertex and edge models with at most 4— body interactions with q' < q. 
Again, only a polynomial overhead is required. 

• We emphasize that the completeness results presented here are based on explicit 
constructions (which cLFG , ctS mentioned, in most cases also computationally 
efficient), which we illustrate with help of a number of examples. For instance, 
we show explicitly how to express the partition function of a 3D Ising model as the 
partition function of a polynomially enlarged 2D Ising model with magnetic fields. 

This paper is organized as follows. First we briefly review the quantum formulations 
of the partition function presented in [21] (section [2]). Then we present completeness 
results with complex parameters (section [3]); that is, we show that the 2D Ising model 
is complete with complex parameters for all models, and that a 2D two-state vertex 
model is also complete with complex parameters for all models. In section [5] we present 
completeness results with real parameters: first we show that the 3D Ising model without 
magnetic fields is complete with real parameters for Ising models on arbitrary graphs, 
then we show that a 3D q— state edge model is complete with real parameters for q'— state 
edge models with pairwise interactions, with q' < q, and finally we prove that a 2D 
q— state vertex model is complete with real parameters for all q'— state edge and vertex 
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models with at most 4— body interactions, with q' < q. In section E] we give some 
implications of our results in complexity theory, and finally we conclude in section [7J 

2. Quantum formulation of the partition function 

Here we briefly review four mappings thoroughly exposed in [21] that allow us to 
express the partition function of a system of classical spins in a quantum language; 
more precisely, as the overlap between a stabilizer state and a complete product state. 
We present mappings for the partition function of the Ising model with magnetic fields 
(introduced in [1]), a general class of q— state edge models (introduced in [7J), both of 
them in the Graph picture, edge models in the GHZ picture and vertex models in the 
PEPS picture (both introduced in [2"T]). 

2.1. Edge models with magnetic fields in the Graph picture 

Generally, an edge model represented on a graph is a model in which particles are sitting 
at the vertices of a graph G = (V, E) and edges represent interactions between these 
particles, thereby resulting in the Graph picture of the edge model. 

Here we consider a classical spin model involving N q— state spins (si, s% . . . s_/v) = s > 
where s a = 0, 1, . . . , q — 1 and a = 1, . . . , N. The spins interact pairwise according to 
an interaction pattern specified by a graph G = (V, E) with inhomogeneous coupling 
strengths J a b and local magnetic fields h a as 

H G (S) = - Jab(25(s a ,S b )-l)-J2haS a , (1) 

{a,b}£E a£V 

where 5(s a , s&) = 1 if s a = and it is otherwise. This corresponds to the Hamilton 
function of the Potts model with magnetic fields, which specializes to the Ising model 
with magnetic fields for q = 2, as considered in [I]. In the remainder of this subsection 
we shall focus on the case q = 2. The partition function Z G is defined as 

z G ({j a6 ,u)=E e ~^ G(s) > ( 2 ) 

s 

where f3 = l/(fcsT), with ks the Boltzmann constant and T the temperature. 

To express this partition function in a quantum language, we first associate with 
each graph G its decorated version by adding a new vertex at the middle of every edge 
(see figure [TJ(b)). We call this new set of vertices "edge qubits" in contrast to the original 
set of vertices, called "vertex qubits" . An edge qubit placed at edge (a, b) e E is labelled 
ab, thus indicating the two vertex qubits a and b which are its neighbors. We now define 
a stabilizer state \<^>g) on the decorated graph: 

l^> = E<8>K + s >><8)K> ( 3 ) 

s ab a 

(for notational transparency we omit the particle outside the ket, that is, we write \s a ) 
instead of \s a ) a ). A stabilizer state is the common unique eigenstate of a set of (Pauli) 
operators; in particular, the state \(pa) is left invariant by the operator ai^ &)b-(a b)&E a ^ 
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for all a G V and by a z a a z a z for all (a, b) & E, where a x and <r^ are the Pauli matrices 
(see [21] for the details). Thus, by applying Hadamard rotations to all edge qubits one 
obtains a graph state (as defined in, e.g., (22J 122]) on the decorated graph. We see 
that this stabilizer state assigns to the vertex qubits the value of the classical spin, 
and to edge qubits the sum of the values of the two neighboring vertex qubits, and 
sums over all spin configurations. In other words, classical spins are associated with 
the vertices of a graph, and we place quantum degrees of freedom at the vertices of the 
graph (vertex qubits) to characterize the degrees of freedom of the local magnetic fields 
of the corresponding classical spins, and quantum degrees of freedom at the edges (edge 
qubits) to characterize the interactions between classical spins. 

We further define a complete product state \ a) = (££), a b ^ eE \a a b) ^S> a ev \ a a)i i n which 

\a a ) = e /3/la |0) +e- ph *\ 1) (4) 

refers to vertex qubit a and contains the information of the magnetic field at that site, 
and 

\a ab ) = e 0Jab \O) + e~ f3Jab \l) (5) 

refers to the edge qubit ab and contains the information of the coupling strength between 
spins at vertices a and b. Here, |0) and |1) are defined as the eigenstates of a z with 
eigenvalues +1 and —1, respectively. 

The partition function of this system of classical spins can be expressed in a 
quantum language by computing the overlap 

Z G = (a\<p G ). (6) 




g Wg) m 

(a) (b) (c) 

Figure 1. (a) A graph G = (V 7 E) is used to indicate that classical spins sit at the 
vertices of this graph and interact with other vertices to which they are connected by 
an edge, (b) The quantum state \(pc) is obtained by placing quantum spins at the 
vertices of this graph, thereby characterizing the local magnetic fields acting on the 
corresponding classical spin, and quantum spins on the edges, thereby characterizing 
the interaction between classical spins along that edge, (c) The quantum state IV'g) 
associated with the graph G places quantum spins only at the edges of G, used to 
characterize the interactions of the classical spins. 



2.2. Edge models without magnetic fields in the Graph picture 

In [7], and in more detail in [21], a similar class of spin models is considered, but there 
the q— state particles are not subjected to magnetic fields, and they interact pairwise 
according to a more general Hamilton function. More precisely, an oriented version of 
a graph G is considered, G a , where, for every edge e, one end- vertex i>+ is assigned to 
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be the head of e, and the other end-vertex v e is the tail of e. The Hamilton function of 
the system is then given by 

H G (s) = Y,h e (\s vt -s v -\ q ). (7) 

eS-B 

Here every h e is some local Hamilton function defined on the edge e with the only 
restriction that the interaction strength between the spins on the endpoints v£ and v~ 
is a function of the difference between s„+ and s„- modulo q, denoted by \s,.+ — s„-L. 
Note that the pairwise interactions here are more general than in ([TJ) , and they can also 
specialize to the Potts model or the Ising model. 

We define a stabilizer state \ip) = Xlsez™ *S>ab I \ s a ~ s b\q)- In order to write its 
stabilizers, we need the definition of the shift and phase operators, X and Z, respectively: 

X\j) = \j + lmodq), and Z\j) = e 2 ^\j), (8) 

for every j = 0,...,q — 1. Then the stabilizer group of \ip G ) is generated by 

a,b)eE -^-ab f° r a ^ a e an d ®(a,b)edf ^ab f° r an f aces / (by face we mean a 
plaquette with the minimal number of edges in its boundary). Here T a b = 1(— 1) if cl 
is the head (tail) of edge ab, (a, b) G df means that edge ab is at the boundary of face 
/, and v a b = 1(— 1) if the oriented edge points clockwise (counterclockwise) according 
to the center of the face (see [21] for a more detailed analysis). We see that the state 
\ipc) is obtained by placing q— dimensional quantum spins at the edges of the graph G, 
which contain the degrees of freedom of the interaction between the two classical spins 
at the end points of the edge (see figure QJc)). 

We also define a complete product state \a) = <S> e eE l a e), where 

KHxy^ij) (9) 

3=0 

is a q— dimensional state associated with the edge e, which contains the information of 
the interaction along that edge. The partition function Zq is obtained by computing 
the overlap 

Z G = q-^ q G \a). (10) 

Relation between \tp) and Note that the state \iPg) specializes to the state \<Pg) 
if the interaction between spins s a and s& connected by e is such that h e (\s v + —s v - \ q ) = 1 
if |s^+ — s v - \ q = 0, and h e (\s v + — s v - \ g ) = otherwise, for all edges e G V. It must also 
hold that the magnetic fields at every site a are zero, i.e. h a = 0, for all a G V . 

However, we can also incorporate magnetic fields in the description of in the 
following way. Consider the state \ip) defined on a certain graph G = (V, E), \i/>g)- Now 
add an extra vertex vq which is connected via decorated edges to all other vertices in G, 
thus forming the state denoted by \ipG+h) (see figure [2]). By letting v interact with all 
other qubits with a coupling strength equal to the magnetic field, i.e. Jo a = h a for all 
a G V, the resulting state is \ipG+h) = 2|v?g) (see [21] for a proof). Note that throughout 
this paper we deal with non— normalized states. Notice also that the number of particles 
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in \(p G ) and \ip G ) is the same: while \ipc) contains \E\ edge qubits and \V\ vertex qubits, 
\ipG+h) contains \E\ + \V\ edge qubits. 



Vq 



\<Pg) \^G+h) 

Figure 2. The state |c/?g) is equivalent to the state \ipG+h)- The latter is obtained 
by taking the state |-0g), adding a new vertex vq connected to all other vertices in G, 
and letting the interaction between vq and every vertex a €z V equal the magnetic field 
at a. 



2.3. Spin models in the GHZ picture 

Here we present the GHZ picture of an edge model [21]. Consider first a q— state edge 
model on a 2D square lattice. We can represent it in its GHZ picture by substituting 
every classical particle sitting at a vertex by 4 (virtual) quantum particles in a GHZ 
state of q— levels: |0000} + |1111} + . . -+\q— 1, q—l, 9—1, q — 1), and by substituting every 
interaction by a projection onto two quantum particles (each belonging to a different 
GHZ state) onto the state \xab) = Z~2 Sa s b e~ /3h( - Sa ' Sb ^\s a Sb) , where h(s a ,Sb) is the part of 
the Hamilton function that describes the interaction between s a and (see figure [3]). 

Then, the partition function of the edge model can be found by computing the 
overlap 

Z e ,2D = (X\GHZ) (11) 

where e indicates that it is an edge model, and 2D indicates that the interaction 
pattern is a 2D square lattice. Here \x) is the tensor product of all the projections, 
Ix) — ®(a6)eE \Xab) (where {a, b} are edges in the graph picture of the edge model), 
and \GHZ) is the tensor product of the GHZ states associated with each particle, 
\GHZ q ) = ® a \GHZ q )&. 

For example, the GHZ picture of this edge model specializes to the Ising model 
without magnetic fields by setting q = 2 and \Xab) = e~ /3Jaf, (|00) + |ll))+e /3Jai, (|01) + |10)) 
for all edges (a, b) in the 2D lattice. Local magnetic fields can be incorporated in this 
picture by setting \ Xa b) = El a ,s b =o exp{/3[(-l) s "+^ J a6 + (-l) Sa h a + (-l) s *h b ]}\s a s b ). 
Similarly, the Potts model on a graph G = (V, E) on q— level systems is obtained by 

setting |xo6> = e^ Jab J2%1 + 2~2j^ m \j m ) for a11 ( a > b ) e E - 

More generally, we can obtain the GHZ picture of a general class of edge models in 
the following way. Consider a q— state edge model in which spins interact according to 
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Figure 3. GHZ picture of a q— state edge model on a 2D square lattice. Rhombi 
represent GHZ states of q— level systems, each of which involves n a (virtual) particles. 
Ellipses represent projections, which describe the interaction between different GHZ 
groups. Here all interactions are pairwise. Notice that the lines between the particles 
in the GHZ state have a different meaning than the edges between vertices in a graph 
state (which correspond to a controlled phase gate between two qubits). 



the Hamilton function 

fr(s) = -5>(s»), (12) 

where % is an interaction (possibly involving many particles) within an index set J, hi is 
the Hamilton function associated with that interaction, and are the spins involved 
in interaction i. 

In order to obtain the GHZ picture of such an edge model [2~T] . we replace each 
particle a by n a (virtual) quantum particles in the state 
9-1 

J2\i)® na = \GHZ q )( a \ (13) 

i=0 

Here n a is the number of interactions in which particle a is involved (e.g., for a graph, n a 
is the number of adjacent edges it has). The state formed by all particles is the product 
state of the GHZ state associated with each particle, viz. 

\GHZ q ) = ($\GHZ q )M. (14) 

a£V 

This is a stabilizer state, and the generators of its stabilizer group are z[ a ^ for 
j E {2, . . . , n a } and <S)j=i -^-j for all a G V, where X and Z are the generalized Pauli 
operators of (jSJ). 

We also replace each interaction i by a projection onto the state 

| Xi) = ^ e ^(s«)| sW)) (15 ) 

which involves |{s^}| quantum particles, each of which belongs to a different GHZ 
state (see figure H] for an example of how to depict many body interactions in the GHZ 
picture). Note that the interaction pattern is encoded in these projections, since they 
specify which spin interacts with which. The dimension d of the interaction particle 
corresponds to the number of different energy values that are assigned. Thus, for 
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A;— body interactions between q level systems, the most general projection has dimension 
d = q k . Note that such projection includes all possible interactions between 2, . . . , k — 1 
bodies, as well as local terms. We define \x) as the tensor product of all projections 
Ix) = ® tg/ \Xi)- Then the partition function of such an edge model in the GHZ picture 
is found by computing 

Z eiI =( X \GHZ q ). (16) 




Figure 4. A Hamilton function involving 2, 3 and 4— body interactions in the GHZ 
picture. Projections contain the information of the interaction pattern (which classical 
spin interacts with which, and whether many-body interactions are present), as well as 
the details of the interactions. Not all interactions are represented in the figure, that 
is why some quantum particles (black dots) are not projected. 



2-4- Spin models in the PEPS picture 

Here we present the PEPS picture of a vertex model (PEPS stands for Projected 
Entangled Pairs) [2H EH ESj- A vertex model represented on a graph G = (V,E) 
is a model in which particles (possibly q— dimensional) are sitting at the edges of the 
graph and vertices represent (typically many-body) interactions between particles. This 
is the Graph picture of a vertex model. 

In order to obtain the PEPS picture of a vertex model, we proceed similarly as for 
the GHZ picture of an edge model. Consider first as an example a vertex model on a 
2D square lattice. Its PEPS picture is obtained by substituting every particle e by a 
Bell pair X^=o =: \4>q)^ (i- e - a GHZ state of two particles), and every 4— body 

interaction by a projection \\ v ) onto one particle of each of the four Bell pairs that 'end' 
in that vertex (see figure [5]). Similarly, the state formed by all particles is the product 
state of the Bell state associated with each particle, viz. |</>+) = <S> eeE \(pg)^ e \ and we 
define \x) as the product state of all interactions, \x) = (£) v ev \Xv) Thus, the partition 
function of this vertex model in the PEPS picture is obtained by computing 

Z vJ =(x\rt)- (17) 

For example, this vertex model can specialize to widely studied models such as 
the six-vertex model or the eight-vertex model. To do this, we express the 4 indices of 
a 4— body interaction as the 2 rows and 2 columns of a matrix, and we write \x v ) as 
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\Xv) 

Figure 5. PEPS picture of a vertex model on a 2D square lattice. 

Xv = ^2ijki e Si ' Sj ' Sk Then this model specializes to the eight-vertex model 

by setting the function h(si, Sj, s^, s/) such that 

/ e? El 
e^ 2 eP E * 
Xv ~ e^ 6 eP* 
K e^ Ea 

and the six- vertex model is obtained by additionally setting e^ Er = e f3Ea = [26J. 
This general class of vertex models can also specialize to the topological colour codes 
introduced in [6], if we let the spins interact according to an hexagonal lattice, and we 
set the Bell states to have = 2 and the interactions \y v ) = Y] a „ „ e ( -~ 1 ^ a+Sb+! " : ^ J \s a + 
Sft + s c ) for all vertices v in the lattice. 

More generally, a PEPS picture can be seen as a specific case of a GHZ picture. 
That is, a picture of a spin model in which particles are substituted by GHZ states, 
and interactions by projections onto particles of different GHZ states. In the previous 
subsection we have shown how to obtain the GHZ picture of a general class of edge 
models. With the same procedure we would describe a general class of vertex models, 
that is, vertex models on q— level systems and with any combination of 1, 2, . . . , k many- 
body interactions, for any q and k, in the PEPS picture. 

3. Completeness of the 2D Ising model with complex parameters 

In this section we prove the first main result of this paper, that is, that the 2D Ising 
model is complete, i.e. that the partition function of any model of classical spins with 
any interaction pattern can be expressed as a special instance of the partition function 
of the classical Ising model with magnetic fields on an (enlarged) 2D square lattice. 
To prove it, we first show that it is complete for Ising models with magnetic fields on 
arbitrary graphs, then for edge models on q— level systems with k— body interactions (for 
any q and k) and finally for vertex models on q— level systems with k— body interactions 
(for any q and k as well). We also show that the enlargement of the 2D square lattice 
is polynomial for Ising models on all graphs, and for edge and vertex models with q 
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and k bounded. However, in order to obtain these results, one needs to allow for some 
coupling strengths and magnetic fields on the 2D Ising model to be complex, and thus 
not physical. 

Finally, we will present another completeness result with complex parameters, 
namely that the 2D vertex model with general interactions is complete with complex 
parameters for all edge and vertex models (including Ising models on arbitrary graphs). 

3.1. Ising models with magnetic fields 

To prove the completeness of the 2D Ising model for Ising models on arbitrary graphs, 
we make use of and of the knowledge of one-way quantum computation. The one- 
way quantum computer is a paradigm for quantum computation, where computations 
are realized by performing single-qubit measurements on a highly entangled substrate 
state called the 2D cluster state \C) [?]; the latter is a graph state [22] associated with 
a 2D square lattice C. 

In particular, we use the fact that the one-way quantum computer is universal. 
This means that any n— qubit quantum state can be prepared, up to local unitary Pauli 
operations, by performing sequences of single-qubit measurements on a d x d cluster 
state \C) of sufficiently large system size M = d 2 . This implies that every n— qubit 
quantum state can be written in the following way: 

E|^) = 2< M - n >/ 2 (/ ® {(3\)\C). (19) 

This formula represents one measurement branch of a one-way computation performed 
on an M— qubit cluster state, yielding the state \ip) (up to a local operation E) as an 
output state on the subset of n qubits which have not been measured. The product state 
(P\ = om y ac ts on the measured qubits, and is determined by the bases and 

the outcomes of the different steps in the computation. The local unitary operator E, 
called "correction operator", acts on the unmeasured qubits (i.e. on the Hilbert space of 
|"0))- The tensor factors of E are always instances of Pauli operators, E, £ {I, X, Y, Z}. 
The prefactor 2 {M ~ n)/2 reflects the fact that the success probability of every measurement 
branch is 2 n ~ M . 

The state \<P2d) (i- e - the state \cpa) where G is the 2D square lattice) is 
also a universal resource. This is because the 2D cluster state \C) can be 
deterministically generated from |<^2d) (up to a local correction) by performing single- 
qubit Y — measurements on all edge qubits. As a consequence, one has 

E / |C) = 2l E l/ 2 (/®(0y| y -)|^), (20) 

where |0y) yj3 is a tensor product of the (+1)— eigenstate of Y on all edge qubits, and 
E' is a local correction. 

We are now ready to establish the connection between the evaluation of Ising 
partition functions and universal MQC To do so, consider the following procedure. 
First, the partition function of an Ising model on a graph G can be expressed in the 
form (jfjj). The stabilizer state \<P2d) is obtained from the cluster state \C) after applying 
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a certain measurement pattern \/3) (i.e. we consider fU9p with = |<^2d))- Finally the 
cluster state \C) is obtained from the state |y?2r>) after measuring all edge qubits in the 
Y basis (see equation (1201) and figure E]). This means that the partition function of the 
Ising model on the graph G can be written as 

Z G ({Ja b ,K})=A( 1 \ip 2D ), (21) 

where A = 2^ E \ + \ v \ +M - n ^/ 2 is a constant and (7) is a product state, (7) = S|a) ®T!\(3) <g> 
lOy)^. Now, by comparing the right hand side of (12T1) with ([6]) we see that it corresponds 
to the partition function of the Ising model on a 2D square lattice evaluated with a set 
of parameters {JL,^} determined by I7). This allows us to conclude that Zq can be 
written as follows: 

Z G ({J ab ,h a })ocZ 2B ({J' ij ,h' i }). (22) 

In other words, the partition function of the Ising model on an arbitrary graph can 
be recovered as a special instance of the partition function of the Ising model on a 2D 
square lattice. This proves that the 2D Ising model is complete for Ising models on any 
other graph. We now give a few remarks regarding this construction. 
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\<Pg) \C) \<P2d) 

Figure 6. The graph state \<pg) (a specific instance of it is drawn in the figure) is 
generated from the cluster state \C) by applying a measurement pattern (/3| on the 
cluster. In its turn, the cluster state is generated from the state |</?2d) by measuring 
all edge qubits Ve (denoted by red dots) in the Y — basis. Black dots denote vertex 
qubits. 



Interaction pattern and details of the interaction. Note that in (|2ip the product 
state I7) is determined by both the interaction graph G and the couplings {J a b,h a } 
of the original model. For, on the one hand, it contains the states \a a b) an d \ot a ) 
encoding the couplings of the original model; on the other hand, I7) contains states 
\j3j) and |0y) corresponding to the sequence of one-qubit measurements which are to 
be implemented in order to generate \^>g) from the universal resource \^>2d)- In going 
from (f2~T|) to (122|) . the state I7) in turn determines the couplings with which the 2D 
model is to be evaluated. Note that the decorated cluster state |<^2d) has vertex qubits 
and edge qubits. The factors of (7) acting on the edge qubits determine the pairwise 
interactions J-j, whereas the factors of (7) acting on the vertex qubits determine the 
external fields h[. The tensor factors of I7) which act on the edge qubits are all equal 
to |0y) oc |0) + i|l). This implies, in particular that in (|22|) only homogeneous pairwise 
couplings JL need to be considered. Hence, all information regarding the couplings J a b, 
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the external fields h a of the original model and the graph G of this model will be encoded 
in the factors of I7) acting on the vertex qubits, and thus in the external fields h\ (which 
will typically be inhomogeneous) . This fact does no longer hold in our construction for 
the completeness results with real parameters (section [5]). 

Polynomial enlargement of the system size. We show that the number of qubits 
in the 2D square lattice is polynomially enlarged w.r.t the number of qubits in the final 
graph. As proved in [19], for all n— qubit states that can be efficiently prepared in 
the circuit model, i.e. by a polynomial sequence of two-qubit gates, the required size M 
of the cluster state in (fl9l) scales polynomially with the number of qubits: M oc poly(n). 
Moreover, in this case the measurement bases \f3j) as well as the correction operations 
E can be efficiently determined. Since any graph state on n qubits can be prepared 
using at most 0(n 2 ) controlled-phase gates [22] , it follows that an arbitrary n— qubit 
graph state can be written in the form ( fT9l) with M = poly(n). Furthermore, for the 
preparation of graph states every single-qubit state \[3j) can always be chosen to be one 
of the X—, Y — and Z— eigenstates. In conclusion, since Ypg) is a stabilizer state, the 
system size of the 2D cluster state grows polynomially with the size of G. Furthermore, 
1/3) consists of X—, Y — and Z— eigenstates. 

Universality of triangular, hexagonal and Kagome lattices. Note that, in the above 
sequence of arguments, one step is particularly crucial, the universality of the 2D cluster 
states: this property is used to "map" an arbitrary state \<Pg), an d hence the associated 
partition function, to the 2D cluster state, i.e all states can be "reduced" to this single 
structure. It was proven in [271 EE] that not only the cluster state is universal, but graph 
states on a triangular, hexagonal and Kagome lattices are also universal resources. 
Therefore, the conclusions for the 2D square lattice apply also to any of these three 
lattices, i.e. the Ising model on a 2D triangular, hexagonal or Kagome lattice is complete 
as well. Also a 2D square lattice with holes has been proven to be universal, and hence 
the same completeness results hold for it as well [29|. 

Complex couplings and fields. In order to interpret measurements with couplings 
(or fields), one needs to identify the coefficients of the projection when expressed in the 
eigenbasis {|0), |1)} with the coefficients e /3J ' , e _/3J (or e^ h , e - ^), respectively. In this 
manner we unambiguously associate a projection with a coupling strength (or magnetic 
field). 

However, the coefficients e ±/3J are real and positive for all (real) values of J 
(similarly for h). This means that if the coefficient of the projection is negative or 
complex, this will translate to a complex value of the coupling strength J or the magnetic 
field h. For example, a projection of an edge qubit onto the eigenstate of Y with 
eigenvalue +1, \0y) oc |0) + i\l) corresponds to a coupling strength J-j = — in/ (4(3) (or 
the same value of the magnetic field if applied on a vertex qubit). Note that we have 
used Y measurements in the proof above in order to transform the state \(£2d) into \^>g)i 
for a general G. More precisely, we have applied Y measurements on every edge qubit 
to transform |(/?2d) to the cluster state \C) (thus corresponding to a uniform, complex 
coupling strength between all vertex qubits, i.e. </• • = — m/(4f3) for all edges (i,j) in 
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the 2D lattice), and the measurements on the cluster state to transform it to \^>g) will 
generally involve Y measurements as well (i.e. some magnetic fields will take complex 
values h! i = —m/{&(3)). This means that the 2D Ising model is complete with complex 
parameters for Ising models on any other graph 

Note that even though such complex coupling strengths or magnetic fields do not 
correspond to physical models, considering the partition function as a function with 
complex arguments is commonly done, e.g., in the context of evaluating the Tutte 
polynomial or finding (complex) zeros of Zq to identify phase transition points [30] . 

A simple analysis reveals when real parameters arise: within Pauli measurements, 
the projectors |0), |1) and |+) = |0) + |1) correspond to real parameters, i.e. X and 
Z measurements (recall that we only consider one branch of the outcomes); and in the 
Bloch sphere, projections onto cos(#/2)|0) + sin(#/2)|l), with 9 G [0,7r] correspond to 
real parameters. However, tilted measurements performed on a stabilizer state typically 
do not yield another stabilizer state [22], and our target state is a stabilizer state. 
Therefore, we shall only consider Pauli measurements. This will be important in the 
next sections, where we address the question what completeness results can be obtained 
with real parameters. 

3.2. All edge models 

Here we show that the partition function of the Ising model on a 2D square lattice can 
specialize to the partition function of a model in which the interacting particles are 
q— dimensional systems. Moreover, we prove that the 2D square lattice is polynomially 
enlarged w.r.t the graph on which the q— dimensional model is defined as long as q 
is bounded. Even more generally, we allow for A;— body interactions and we find that 
these models can also be generated efficiently provided k and q are bounded. Even 
though an unbounded growth of q or k with the system size may require an exponential 
enlargement of the system size, we are able to give some examples in which either q or 
k is not bounded and we still obtain a polynomial enlargement of the system size. 

In order to prove these results, we use the GHZ picture presented in section 12.31 
We want to express the partition function of ( fl6l) as a special instance of the partition 
function of the 2D Ising model. To do so, we observe that the projection onto the state 
\Xi) can be implemented by first preparing the state \x%) on the quantum particles on 
which the projection is to be performed (i.e. qubits belonging to the different GHZ 
states), and then measuring locally these quantum particles in the X basis. The 
preparation of the state \xi) (of dimension d) can be done by a unitary gate U\ acting 
on = \log 2 d^] qubits, \xi) = Ui\0)® md (see figure [7] for a two— body interaction, 
and figure [8] for a 4— body interaction)). Thus, the product state of all interactions can 
be expressed as |%) = U\0)® md where U = &) ieI Ui and m d = ^ie/ m i ■ This implies 
that the partition function of a q— state model can be written as 

Z q G = (GHZ q \U\0)® md . (23) 

The state |0)® md is defined on qubits, and the state (GHZ q \U can be prepared from 
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a cluster state by applying a certain measurement pattern on it. This proves the 
completeness of the 2D Ising model for q— state models. Note that the preparation 
of the state (GHZ q \U will generally require measurements that contain the complex 
unit "i" or negative values, and thus the associated coupling strengths or magnetic 
fields in the 2D square lattice will be complex. 




Figure 7. A projection \xi) can be seen as a gate Ui that prepares this state on the 
qubits involved in the interactions, and then these qubits are measured locally in the 
X— basis. For a projection of dimension d, the gate acts on = |~log 2 d] qubits, 
2mj^ = 2[log 2 d] qubits in the one-way model, q— dimensional particles are viewed 
as m q two-dimensional particles. The gate Ui can be decomposed into one— and 
two— qubit gates, which correspond to a certain measurement pattern on the cluster. 
The dotted lines indicate variable size. 

Efficiency results. We now study whether the 2D square lattice on which we 
compute the partition function has to be polynomially or exponentially enlarged w.r.t 
the graph on which the q— state model is defined. 

We first observe that the state \GHZ q ) can be efficiently prepared from a cluster 
state because it is a stabilizer state [22]. In fact, every \GHZ q )^ can be formed of 
m q = |~log 2 q] \GHZ 2 ) states, simply by regarding the values of the latter as the binary 
expression of the q levels that \GHZ q ) can take (see figure Cj). Taking into account that 
the state IGHZq)^ contains n a particles, the preparation of this state requires 0(n a m q ) 
qubits. Also the preparation of a GHZ state of n particles requires n — 1 phase gates. 
Hence the number of gates required to prepare \GHZ q )^ scales as 0(n a m q ). 

We now consider the efficiency in the generation of the unitary gates Ui. Recall that 
each of these gates acts on = [log 2 d^] qubits initialized in |0), where d is at most 
q k for k— body interactions between q— level systems, i.e. < k\og 2 q. These gates can 
be implemented in the model of one-way quantum computation [19] (see figures [7] and 
E]). Note also that whenever two projections cross each other, in the Graph picture this 
corresponds to implementating the gates specific to each projection and then applying 
several SWAP gates that exchange the position of the upper and lower lines and thus 
produce a crossing (see figure [9]). Any such gate can be decomposed into single and two 
qubit gates, each of which requires a constant number of qubits to be implemented in 
the one-way model. Hence we need to study the number of single and two qubit gates 
required in the decomposition of \Xi)- m general an exponential number of gates might 
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be required to prepare the state \xi)- in the following we study which states \x%) require 
a polynomial number of gates. 
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Figure 8. A projection (of dimension d) can be seen as a gate Ui acting on 



Xi) 



= [log 2 d] qubits. g— level particles are viewed as m q two— level particles. 




Figure 9. A crossing of projections in the GHZ picture (right figure) is seen in 
the Graph picture as a the application of the gates corresponding to each projection 
followed by a SWAP gates that produce the crossing (left figure). Auxiliary qubits are 
not represented in the gates for the sake of clarity. 



Let us first assume that q is bounded. A Hamilton function describing k— body 
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interactions contains at most 

(24) 

terms. We have already noted that, if this k— body interaction has its maximum 
dimension d = q k , this includes all possible 1, 2, . . . , k — 1 interactions between these k 
particles. If we further require that k is bounded, then the unitary Ui can be generated 
from a poly-sized circuit, and thus it requires a polynomial number of gates in the one- 
way model. It follows that any model with q and k bounded can be prepared efficiently 

If k is unbounded, the number of terms in the Hamilton function may grow 
exponentially, and the gate t/j may be acting on an unbounded number of qubits. On 
the other hand, if q is unbounded the gate Ui may be acting on an unbounded number 
of qubits as well. Therefore, models with k or q unbounded will generally require an 
exponential overhead in the system size. 

We can find some exceptions to this result. In particular, we show that models 
in which remains bounded even though q or k grow unboundedly can still be 
generated efficiently. More precisely, we want to see that the Potts model with pairwise 
interactions [1] but with q unbounded can still be generated efficiently The idea is that 
the Potts model only assigns two values of the energy ( J e if they are in the same spin 
state, and otherwise) independently of the number of spin states available q, and thus 
m q . The projection corresponding to this model is 

|a e ) = (e^-l)(|0)|0) + |l)|l) + ...+ 

|g-l)|g-l»+ \x)\v) (25) 

x,y=0 

for each edge e G E. Each q— level system is formed of m q = |~log 2 q] qubits, which we 
label xi . . . x mq for one system and y\. . . y mq for the other. Then the state \a e ) can be 
written as 



l«e) = (e^ J -l)(g)|^)^ + | + +)^, (26) 
i=i 

where \(pt)x iyi is a Bell pair between particles Xi and yi. Here |+) and |— ) are the 
eigenstates of a x with eigenvalues +1 and —1, respectively. By ordering the particles as 
x\y\X2y2 ■ ■ ■ %m q ym q , one can see that the maximum rank of the reduced density matrix 
of all bipartitions between the first k particles and the rest of the system is 3. This 
means that, in the Matrix Product State (MPS) representation of this state, all tensors 
have dimension 3 independent of the number of qubits m q [31] . Since this representation 
is efficient, the state can be efficiently generated |32j. 



3.3. All vertex models 

We want to show that the partition function of the 2D Ising model can also specialize to 
the partition function of vertex models with q— dimensional particles with an arbitrary 
interacion pattern. 
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In section we presented the PEPS picture of a vertex model, and we noted that 
it is a particular case of the GHZ picture. It follows that the completeness proof is 
completely analogous as for edge models in the GHZ picture. We conclude that the 
partition function of a vertex model on q— level systems can be written as 

Z q G>v = «l^|O) 0md , (27) 

from where the state jO)®" 7 ^ is defined on qubits, and the state (<f>t\U can be prepared 
from a cluster state by applying a certain measurement pattern on it (see figure [10] for 
a 4— body interaction). This proves the completeness of the 2D Ising model for vertex 
models on q— level systems, for arbitrary q and with an arbitrary interaction pattern. 
Note that the preparation of the state (<f>q\U will also generally require measurements 
that contain the complex unit "i" or negative values, and thus the associated parameters 
will lie in the complex regime. 




Figure 10. A projection onto a state \x%) can be seen as a preparation of this state 
by the gate Ui and the measurement of the particles in the X— basis. 

Efficiency results. Since the PEPS picture is a special case of the GHZ picture, it 
follows that the efficiency results are the same. That is, whenever q and k are bounded, 
the 2D square lattice requires a polynomial enlargement w.r.t to the graph where the 
vertex model is defined. A model with q or k unbounded may require in general an 
exponential overhead. 

3.4- Completeness of the 2D vertex model with complex parameters 

Here we want to show that a vertex model with two— level systems (q = 2) defined on 
a 2D square lattice and with one- dimensional four body interactions is complete with 
complex parameters for the same models for which the 2D Ising model is complete. 



CONTENTS 



20 



The proof is based on the following observation: the PEPS picture of a vertex 
model on a 2D square lattice with two— level systems with one— dimensional projections 
(as the one depicted in figure [5]) corresponds indeed to the PEPS picture of a 2D 
cluster state [2U [25] . More precisely, a single qubit measurement in a 2D cluster state 
would correspond in the PEPS picture to one-dimensional projections with a Hadamard 
rotation on two of the four qubits, since the edges in a cluster state correspond to the 
states |0+) + 11— ), whereas the pairs (which now act as edges) are in the state 1 00) + |11) 
(note that |±) oc |0) ± |1)). This means that all states that can be prepared from a 
cluster state (with measurements corresponding to complex parameters) can also be 
prepared from a 2D two— state vertex model (with complex parameters as well). Since 
the 2D cluster state is universal, it follows that a 2D two— state vertex model is complete 
with complex parameters for all q— state edge and vertex models. Note that all efficiency 
results for the 2D Ising model also hold for this 2D two— state vertex model. 

4. Example: Partition function of the 3D Ising model obtained from the 
partition function of the 2D Ising model 

As an explicit example of the construction presented above, we show how to compute 
the partition of the Ising model with magnetic fields on a 3D square lattice (3D Ising 
model) by computing the partition function of the Ising model on a (polynomially) 
enlarged 2D square lattice. To do so, we only need to find a measurement pattern 
(P\ such that, when applied on the state \<P2d), the resulting state is |<^3d), i-e. the 
stabilizer state \ip) defined on a 3D square lattice. In principle we could generate |</?3d) 
by implementing gates (consisting of measurements) on the cluster state. We present 
here a more direct procedure that requires a smaller cluster to generate the same \<P3d), 
and that is translationally invariant (if the initial resource state is \tp) on a 2D square 
lattice with crossings, instead of |v?2d))- 

This procedure is based on measurement rules, which translate into a certain 
measurement pattern (7] applied on the initial state defined on a certain graph G (e.g. 
\^Pg))i transforming it to the same state defined on another graph G' (e.g. \<Pg'), with 
\<Pg')) = (I <8> (7I) I </?(;}). It is clear from this definition that the measurement rules 
depend on the state they act upon. 

Here we will make use of the following measurement rules for the state \ip) (see ([3]) 
for the definition of | </?)): 



where G' is the same graph as G but with the edge a'b' merged (or: "contracted"), G" is 
the same as G but with the edge a'b' deleted, and G'" is the same as G but in which two 
pairs of particles now have a crossed interaction and other particles have been deleted 
as specified below. 



((+a'|(0 a '6'| ® I)\<Pg) = \<PG') (merge rule) 
{{+a'b'\ <8> I)\<Pg) = \<PG") (deletion rule) 
((7I ® I)\<p G ) = \<Pg»>) (crossing) 



(28) 
(29) 
(30) 
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A merge rule is a measurement pattern according to which two vertex qubits and 
the edge qubit in between are merged into one vertex qubit. The neighbors of the 
resulting vertex qubit are the neighbors of the two vertex qubits that have been merged 
(see figure EDJa)). Several merge rules can be concatenated, thus merging several vertex 
qubits and the edge qubits inbetween into one vertex qubit (see figure fT2l) . According 
to section I3.1[ we can associate each projection with a value of the coupling strength 
or magnetic field, see equations (jlj) and ([5]). In this way we see that the projection of 
the vertex qubit a' onto (+| corresponds to setting the magnetic field h a i =0, and the 
projection of the edge qubit a'V onto (0| corresponds to setting the coupling strength 
J a >b' = oo. That is, an infinite coupling strength between two particles effectively 
corresponds to having one particle. Interestingly, we also see that this rule only requires 
real parameters to be realized. Whether a rule requires real or complex parameters will 
be a relevant issue for the completeness results with real parameters that we will present 
in section [51 

We would like to make a side remark here. We have just seen that one requires 
to set J a iy = oo in the above-presented rule. One may notice that, when this limit is 
taken, the state e^ Ja ' b '\0) + e~^ Ja ' b ' |1) « e f3Ja ' b '\0) instead of the state |0) is obtained, 
which may seem problematic because of the diverging normalizing factor. Here we note 
that there is a formulation of the present results in which this problem does not occur. 
This involves a rescaled version of the Ising model, considering an interaction which 
assigns a value J a b to non-aligned spins and the value to aligned spins. While being 
equivalent to the formulation of the Ising Hamilton function (TjQ), in this formulation the 
states ([5]) become |0) + e~ /3Jab \l). Hence the state |0) is also obtained by taking the limit 
J a b — > oo, but without the occurrence of a very large constant. It is easy to verify that 
all results in the present paper carry through completely in this rescaled formulation. 
Thus we shall henceforth stick to the conventional formulation of the Ising Hamilton 
function ([[]). 

A deletion rule is a measurement pattern according to which two vertex qubits 
become disconnected (see figure fTTlfb)). A projection of an edge qubit a'b' onto 
(+| = (0| + (1| corresponds to a coupling strength J a i h ' = 0, that is, qubits a' and 
b' are not coupled to each other, which corresponds to having no edge between a and b. 
We see, thus, that a deletion rule is realized with real parameters as well. 

The creation of a crossed interaction between two pairs of particles corresponds 
to creating a crossing of edges. This requires the measurement pattern (7I that is 
specified in figure[13l which involves ^-measurements, and hence corresponds to complex 
parameters. 

In the following we shall prove these three rules. Our strategy will be to consider 
the measurement rules known for graph states (see, e.g., [22]) and translate them to 
the state \<p) (in [21] the modification of the state is studied by looking at the stabilizer 
state). We recall that the state \<pa) can be obtained from the graph state defined on 
the decorated version of G, \G), by applying a Hadamard rotation on all edge qubits, 
Wa) = ® ah H ab \G). 
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(b) 



(d) 



1 



Figure 11. Measurement rules for the state \ip). (a) Merge rule. We denote the three 
qubits that are merged with a thick grey line on them. First, the X measurement on 
the vertex qubit is performed, and then the Z measurement on the chosen neighbor. 

(b) Deletion rule. An edge qubit is deleted by applying an X measurement on it. 

(c) Deletion of a vertex qubit by a Z measurement, (d) Y— measurement on an 
edge qubit and on a vertex qubit. Only rules (a) and (b) result in a decorated 
graph. Measurements (a), (b) and (c) correspond to real parameters, measurements 

(d) correspond to complex parameters. 



I 



Figure 12. Concatenation of merge rules. A thick gray line over certain vertex qubits 
(and edge qubits inbetween) symbolize that they have to be merged. 



The reader may skip these proofs and proceed directly with the example. 

Proof of the merge rule. The merge rule for graph states reads: first X— measure the 
vertex qubit a and then Z— measure the chosen neighbor bo [22]. The X measurement 
requires the correction operator 



<iY 



% 

71 



i 



(31) 



on the chosen neighbor. The chosen neighbor is the edge qubit on which the Z 
measurement is to be performed. This means that a projection of the chosen neighbor 
onto (0| actually becomes (0|(\/?F)t oc (+|, i.e. a projection onto an eigenstate of X. 
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Figure 13. Measurement pattern to obtain a plaquette with a decorated crossing 
starting from a 2D square lattice. It involves Y measurements, and hence it corresponds 
to complex parameters. The figure shows explicitly how the underlying graph changes 
with the measurements applied. 

The X measurement also requires the correction operators Z on the neighbors of a 
which are neither 6 nor neighbors of bo . 

We now translate the merge rule for graph states to a merge rule for the state \ip). 
The X measurement on the vertex qubit remains the same for \tp). The projection of 
the chosen neighbor onto (+| becomes (0|, i.e. a Z measurement, after the Hadamard 
rotation. The correction operators Z on the other neighbors of a become X operators. 

Note that the merge rule (the measurements and the corrections) all correspond to 
real parameters. More precisely, first we regarded projections onto (0|, (1| and (0| + (1| 
as corresponding to real parameters. Now we see that, when taking the corrections into 
account, only the projections (0| and (0| + (1| correspond to real parameters. 

Proof of the deletion rule. In the state \ip) the interaction between two particles a 
and b (the original classical spins) is described by the edge particle ab. Hence, a deletion 
rule corresponds to deleting the edge particle and also the edges between particles a and 
b. 

To delete a particle in a graph state one has to measure its spin in the Z— direction 
[Z— measure it). This measurement (more precisely, a projection on (0|) requires no 
correction operators on its neighboring qubits. Thus the Z measurement becomes an X 
measurement on an edge qubit in the state \tp). 

Proof of the creation of a crossing. In order to show how to create a crossing, 
we need to introduce two measurement rules: the deletion of a vertex qubit, and 
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Y —measurements on edge and vertex particles. The deletion of a vertex qubit consists 
of measuring a vertex qubit in such a way that this qubit becomes disconnected from 
all its neighboring edge qubits. As pointed out above, the deletion rule in a graph state 
consists of a Z measurement, and it requires no correction. Hence, a Z measurement on 
a vertex qubit of \<p) corresponds to deleting that qubit (see figure fTTTc)). Regarding Y 
measurements, for graph states it is known that a Y measurement applied on a particle 
transforms the graph state by applying the local complementation rule to that particle 
(i.e. deleting or adding edges between the neighbors of that particle if the edges were or 
were not there, respectively; see [22]) (see figure EDJd)). 

We combine these two rules in the measurement pattern of figure [13] in order to 
create a decorated crossing (that is, the edges that cross each other are also decorated). 
This concludes the proofs of the measurement rules. 

Now we can proceed with the example. The outline of the procedure is as follows: 
first, we generate a 2D square lattice with crossings, with all edges decorated, from 
a 2D decorated square lattice. This involves creation of crossings, which correspond 
to complex parameters. Then, we generate a decorated 3D square lattice from the 
decorated 2D square lattice with crossings. This involves only merge and deletion rules, 
which correspond to real parameters (see figure [TJ 
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Figure 14. We will first obtain a decorated 2D square lattice with crossings from 
a decorated 2D square lattice. This procedure involves complex parameters. From 
the latter lattice, we will obtain a decorated 3D square lattice, in a procedure which 
involves only real parameters. 



In order to generate a decorated 2D square lattice with crossings from a decorated 
2D square lattice, we need to generate plaquettes with decorated crossings, and 
plaquettes without crossings. To obtain a plaquette with a crossing we proceed as 
in figure [13j To obtain a plaquette without a crossing, we only need to select a square 
in the 2D square lattice of the same size as in figure [T3], and delete all vertices inside. 
Then we merge the vertices at the boundaries so that only one decoration remains at 
each side of the square (see figure fT5l) . 

Now we want to generate a 3D square lattice starting from a 2D square lattice 
with crossings by means of the merge and deletion rule alone. To do so, we first embed 
the figure shown in figure fT6T a). which we call a "face", on the 2D square lattice with 
crossings (the face can be seen as part of a three-dimensional structure, as will be made 
explicit later in figure [T9l) . We do it by tilting every square to the left so that the former 
vertical lines of the squares now coincide with the diagonal lines (going from the upper 
left corner to the lower right corner) of the 2D lattice with crossings (coloured tilted 
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Figure 15. Measurement pattern to generate a plaquette without a crossing. The 
plaquette needs to have the same size as the plaquette with a crossing. This pattern 
consists of deleting all inner vertices and merging the edges at the boundaries so that 
only one decoration remains at each side. 

squares in figure [ToTb)). Only the lines crossing the squares remain to be embedded, and 
these are identified with the diagonal lines going from the lower left corner to the upper 
right corner in the 2D square lattice with crossings (thick black lines in figure [ToT b)). 
As shown in the figure, some vertices have to be merged, and also the sides have to be 
merged so that only one decoration remains at each side. Finally we apply the deletion 
rule to all edges that do not correspond to any edge of the projected 3D square lattice. 

Note that the embedding of larger faces can be straightforwardly derived from this 
one. For example, an addition of squares in the horizontal direction in figure [ToT a) 
would translate to an addition of adjacent tilted squares in (b). Likewise, an addition 
of squares in the vertical direction would translate to an addition of "layers of tilted 
squares" in (b). 

Now we embed the concatenation of faces shown in figure fTTT a) in the 2D square 
lattice with crossings. To do so, we shift the position of one embedded face down and to 
the left w.r.t the other one (figure fTTT b)). In this way the lines that join one face with the 
other naturally follow their way while crossing the squares appropriately. Note that this 
embedding is translationally invariant if periodic boundary conditions are considered. 

Next we consider a figure consisting of 4 concatenated faces, each of them with 9 
squares (figure [T8l(a)) with open boundary conditions. We embed this figure in the 2D 
lattice with crossings as shown in figure fTBlfb). 

It only remains to be shown that such concatenation of faces corresponds to a 3D 
square lattice projected in a plane. It indeed corresponds to a 3D square lattice if one 
projects this lattice so that the faces parallel to the plane of the paper do not cross 
each other, and hence the connections between the faces involve crossings of edges (see 
figure fT9l) . In this manner we obtain a systematic way to identify a 3D square lattice 
with a non-planar 2D figure. 

Thus, we have shown an explicit, constructive way to write the Hamilton function 
of the Ising model on a 3D square lattice by writing another Hamilton function on an 
(enlarged) 2D square lattice with complex parameters, or on a 2D square lattice with 
crossings with real parameters. 

Polynomial overhead in the system size. We want to see that the 2D square lattice 
with crossings is only polynomially enlarged w.r.t the 3D square lattice. To see this, first 
consider a face of squares of dimensions ixz, with the shown in figure [19j Here 
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Figure 16. (a) Face of 4 squares, (b) Embedding of this face on a 2D lattice with 
crossings. A square of a certain colour in (a) corresponds to the square of the same 
colour in (b). All crossings are indicated with an empty circle. Most vertices and 
crossings are labelled with a letter to indicate the correspondence of points between 
(a) and (b). All edges in (a) and (b) are decorated, even though this is not represented. 
Wide gray lines indicate vertices that have to be merged. The edges of the underlying 
2D lattice with crossings that have no other line on top of them have to be deleted 
by deleting its decoration. Finally, edges a to d have to be merged, so that only one 
decoration remains between a to e, and analogously for other sides of the squares. 



x (y, z) is the number of vertices that the 3D square lattice has in the x (y, z) direction 
(not the number of squares). Each of the squares has x — 1 crossings in each of its 
vertical sides (e.g. figure [171(a) shows two faces of squares, each face of dimensions 4x3, 
and where each of the squares involves 3 crossings in each of its vertical dimensions). 
The embedding of such face in a 2D square lattice with crossings (figure UTT b)) involves 
2x + l vertices in the horizontal direction and (x+l)(z — l) vertices in the tilted direction 
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(b) 

Figure 17. (a) 2 concatenated faces of 6 squares each, (b) Embedding of these two 
faces in a 2D square lattice with crossings. The remarks of figure [TBI hold for this one 
as well. This embedding is translationally invariant if periodic boundary conditions 
are considered. 

going from the lower right corner to the upper left corner (counting the merging at the 
upper part of the tilted squares). Thus such face requires (2x + l)(x + l)(z — 1) vertices 
in the 2D square lattice with crossings. 

The concatenation of faces only adds the additional vertices on which thick black 
lines go through without having any face at the background (lower right corner and 
upper left corner in figure fT7Tb)). In particular, every concatenation of 2 faces requires 
(2x + l)(x + 1) additional vertices. Therefore the embedding of a 3D square lattice of 
dimensions x x y x z involves (2x + l)(x + l)[(z — l)y + (y — 1)] ~ vertices in the 2D 
square lattice with crossings. Finally, the 2D square lattice with crossings is produced 
from an enlarged 2D square lattice as summarized in Fig [181(c) . We see that the 2D 
lattice is enlarged by a factor of 4 x 4 = 16 w.r.t the 2D lattice with crossings. Hence, 
we require 

16(2x + l)(x + l)[(z - l)y + (y - 1)] ~ 0{x 2 yz) (32) 

vertices to prepare a 3D lattice of size x x y x z. This proves that the number of qubits 
of the 2D lattice is polynomially enlarged w.r.t the number of qubits of the 3D square 
lattice n = xyz. In particular, if one considers a 3D square lattice with x = y = z, 
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Figure 18. (a) 4 concatenated faces, each of them with 9 squares, with open boundary 
conditions, (b) Embedding of this figure in a decorated 2D square lattice with crossings, 
(c) The decorated 2D square lattice with crossings is prepared from a decorated 2D 
square lattice by implementing the measurement pattern specified in the figure. This 
summarizes the procedure. 



then size of the 2D square lattice grows as n 4//3 . For example, to construct a 3D lattice 
of size x = y = z = A we need a 2D square lattice with 10800 vertices, i.e. a lattice of 
104 x 104 qubits. But this is only the number of vertex qubits, to which one has to add 
the number of edge qubits. Since a 2D lattice with crossings of size x' x x' with periodic 
boundary conditions has 3(x' — l)x' edges, and this side has x' = 104, the required 2D 
lattice with crossings has in total 42936 qubits. 

The procedure to project a 3D square lattice on a 2D plane shown in figure [T9l 
provides a polynomial upper bound in the number of crossings needed to represent a 
3D square lattice in a 2D (non-planar) figure. More precisely, for a 3D square lattice 
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Figure 19. (a) 3D square lattice, (b) Projection of this 3D square lattice in a plane. 
Colourful squares correspond to the same squares in figures [T6j [17] and [TH All edges 
are decorated even though the decorations are not represented. 



of dimensions x x y x z the number of crossings c involved in the projection of this 
figure is c = x(x — \){y — \){z — 1) ~ 0{x 2 yz). For example, the 3D square lattice 
of figure [T9f a) has dimensions 4x4x4, and its projection involves c=108 crossings, 
as shown in figure fTWb). Therefore the number of crossings involved in this projection 
also grows polynomially with the 3D square lattice and with the same exponent 4/3. 

5. Completeness results with real parameters 

In this section we restrict the parameters of the complete model to be real, and thereby 
allow for a physical interpretation. We prove that 

• The 2D Ising model with magnetic fields is complete with real parameters for all 
Ising models with magnetic fields on planar graphs. The overhead in the system 
size is polynomial. 

• The 3D Ising model without magnetic fields (or the Ising model without magnetic 
fields on a 2D square lattice with crossings) is complete with real parameters for all 
Ising models with magnetic fields. The overhead in the system size is polynomial. 

Regarding completeness results with for q— state models, we argue that 

• A 3D q— state edge model (or a 2D 2q— state model) is complete with real 
parameters for all q'— state edge models with pairwise interactions, with q' < q. 
The overhead in the system size is polynomial. 

• A 2D q— state vertex model is complete with real parameters for all q'— state edge 
and vertex models with at most 4— body interactions, with q' < q. The overhead 
in the system size is polynomial. 

5.1. Completeness with real parameters of the 3D Ising model without magnetic fields 
for Ising models with magnetic fields 

Completeness with real parameters of the 2D Ising model with magnetic fields for Ising 
models with magnetic fields on planar graphs. It has been argued in section 0] that 
only the merge and deletion rule for the state \cp) correspond to real parameters, and 
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that there is no combination of them that can create a crossing. Consequently, if we 
start from the Ising model on a 2D square lattice, which is a planar graph, it can 
only prepare planar graphs under merge and deletion rule. It follows that the 2D Ising 
model with magnetic fields is complete with real parameters for Ising models on planar 
graphs. Moreover, the 2D square lattice needs to be polynomially enlarged w.r.t the 
graph on which the Ising model is defined, because the final state, \<£g) with G planar, 
is a stabilizer state, and thus it can be prepared efficiently starting from |y?2Z)) [22] . 

Completeness of the 3D Ising model without magnetic fields for Ising models with 
magnetic fields with real parameters. In section [2~2l we have seen that the state can 
incorporate magnetic fields in its description by letting a new vertex vq interact with 
all other vertices with a coupling strength equal to the magnetic field (figure [2]). In 
order to prove the main results of this section, we need to introduce the merge and 
the deletion rule for the state which will be the only two rules corresponding to 
real parameters. They coincide with the measurement rules for the state \ip) (see the 
beginning of section Hj) but without the vertex qubits and generalized to q— level systems 
(see the end of section [2T21 for the relation between and | </?)). 

Merge rule: Z— measurement. The merge rule for the state \ipc) consists of 
measuring the edge qubit ab in such a way that the resulting state \4>g') has a' and 
b' as the same vertex. It coincides with the merge rule for \tp) in the edge qubit: it 
consists of projecting the edge qubit on (0|, i.e. it is one branch of a Z measurement on 
the edge qubit (see figure 1207a)). Note that this projection corresponds to a coupling 
strength ho = oo and h{i) = — oo for ? ^ 0, that is, when s a / — sy = the coupling is 
infinitely strong and they effectively become the same particle. 

Deletion rule: X —measurement. It coincides with the deletion of an edge qubit on 
the state \<p): it consists of projecting the edge qubit ab onto (+ g \, thus corresponding 
to one branch of an X measurement (see figure 1207b)). Note that this projection 
corresponds to a coupling strength h(j) = for all j, that is, the spins a and b have 
become disconnected. 



(a) 



(b) 
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Figure 20. Measurement rules for \tp). (a) The merge rule consists of Z-measuring 
the edge qubit. (b) The deletion rule consists of X-measuring the edge qubit. Both 
measurements correspond to real parameters. 



Now we can proceed with the main results of this section. Here we want to prove 
that the 3D Ising model without magnetic fields (represented by the state |^3d)) is 
complete with real parameters for Ising models with magnetic fields on arbitrary graphs 
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(i.e. \<Pg) with arbitrary G). 

We first observe that from IV^d) we can prepare the state \ip 3 r,'+h) = 2|y? 3jD /) by 
applying only the merge and deletion rules. This state corresponds to the 3D Ising 
model with magnetic fields on a smaller lattice (see figure |2T1) . To see this, we transform 
a certain part inside of the 3D square lattice to a 2D square lattice with crossings as 
shown in section HI figure [T8l This part will correspond to the smaller 3D square lattice, 
3D'. Then we use the remaining part of the 3D square lattice to create one vertex 
(vq) connected to all vertices of 3D' (black dots in figure [TBI) . Finally the 2D square 
lattice with crossings is converted to a 3D square lattice again. Note that this procedure 
only requires merge and deletion rules. Moreover, it shows that there is a polynomial 
overhead in the number of spins between \ip^D) an d \"^3D'+h)- Hence, it only remains to 
be shown that the 3D Ising model with magnetic fields is complete with real parameters 
for Ising models with magnetic fields on arbitrary graphs, which is to what we turn 
next. 
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Figure 21. The 3D Ising model without magnetic fields (IV^d)) can prepare by 
means of the merge and deletion rules the 3D Ising model with magnetic fields 
WsD'+h) = 2|<^3fl'»- 



Note that, in terms of complexity classes, the 3D Ising model without magnetic 
fields is known to be as hard as the 2D Ising with magnetic fields (it is an NP- 
complete problem). However we see that the completeness results with real parameters 
are different for the two models: while the 3D Ising model without magnetic fields is 
complete for Ising models with magnetic fields on arbitrary graphs, the 2D Ising model 
with magnetic fields can only prepare Ising models with magnetic fields on planar graphs. 

Completeness of the 3D Ising model with magnetic fields for Ising models with 
magnetic fields. We have seen in section H] that the state \<pa) has to be defined on a 
non-planar graph in order to prepare \<Pg') defined on other non-planar graphs if only 
the merge and deletion rules can be applied. A natural candidate for a non-planar 
graph is the Ising model with magnetic fields on a 3D square lattice (3D Ising model) 
which indeed is complete with real parameters for Ising models on arbitrary graphs. 

To prove the latter claim, we first observe that a resource is complete if and only 
if (iff) it can prepare a decorated clique (i.e. a fully connected graph with all edges 
being decorated). The forward implication is trivial, for if a resource can prepare any 
decorated graph, a particular instance of such decorated graph is a decorated clique. 
To see the backwards implication, note that, by definition, a decorated clique contains 
all possible (decorated) edges between the vertices. Thus, in order to obtain a target 
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decorated graph of n vertices one only needs to delete the unnecessary edges of the 
decorated clique of n vertices by deleting its decoration. 

Therefore we need to show that a 3D square lattice can prepare a decorated clique 
only by merging and deleting qubits. To construct a decorated clique with n vertex 
qubits, let us first define the quantities m, n 2 and 77,3 satisfying 2(m + n 2 ) — n — 1 and 
n 3 = 2{n— 4). Now we consider a 3D square lattice of size (3ni+2ri 3 ) x (3n 2 +2n 3 ) x (n— 3) 
(axes shown in figure [T9|) and proceed as follows. We take the lowest layer (x — y plane 
with z = 1) and merge all the vertices inside this layer in order to form a hub connected 
to all the vertices in the boundary via decorated edges (see figure l22Ta)). We then delete 
2 out of every 3 edges that connect the hub to the boundary in a successive way, and 
merge vertices at the boundary so that only the vertices connected to the hub remain. 
Notice that the hub has now 2(ni + n 2 ) = n — 1 neighbors in the boundary, and thus 
this vertex already has all its connections to the other vertices. Now we need to connect 
the n — 1 remaining vertices among each other. 




(a) (b) 

Figure 22. Construction of a decorated clique with n = 7 vertices. Procedure for the 
lowest layer of the 3D square lattice, (a) First build a hub with 3(n — 1) neighbors. 
Delete 2 out of every 3 edges connected to the hub, and merge vertices at the boundary 
so that only the vertices connected to the hub remain, (b) Use the remaining lattice 
around the hub to merge vertices to the boundary qubits so that each of them has 
n — 4 new neighbors. 

To do so, we first merge n — 4 vertices of the remaining lattice of this layer to each 
boundary qubit. Then we delete all other edges in this layer until we obtain figure 1227 b). 
Now we label boundary qubits as 1, 2, . . . , n — 1 and the new neighbors to each boundary 
qubit as 1, 2, . . . , n — 4. The idea is that vertex 1 connects to the other boundary vertices 
through vertices 1, and similarly for vertices up to n — 4. To do so, we connect each 
vertex 1 to the qubit above it in the 2nd layer (x — y plane, z — 2). In the 2nd layer we 
connect them using a hub, all whose connections have to be merged. We also merge each 
vertex 1 to its boundary qubit, and qubit 1 of the boundary qubit 1 to its above qubit 
(see figure [23]). We repeat the same procedure for qubits 2 through the 3rd layer, and so 
on up to qubit n — 4 through the n — 3 layer. Finally we delete the edges that connect 
each pair of boundary qubits more than once and we obtain the decorated clique shown 
in figure EH 
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This construction also shows that the 3D square lattice needs to be polynomially 
enlarged w.r.t the target graph. More precisely, it requires (3ni + 2n 3 ) x (3ri2 + 2n 3 ) x 
{n — 3), with 2{n\ + n?) = n — 1 and n 3 = 2{n — 4), i.e. 0(n 3 ) spins for a target graph 
with n spins. For example, to construct a clique with n = 7 vertices as the one shown 
in figure EH one requires a 3D lattice of size 10 x 9 x 4 = 270 vertices. 

This proves that the 3D Ising model with magnetic fields is complete with real 
parameters for Ising models with magnetic fields, and that the overhead in the system 
size is polynomial. Due to the remarks at the beginning of this Section, this also proves 
that the 3D Ising model without magnetic fields is complete with real parameters for 
Ising models with magnetic fields, and the overhead in the system size is also polynomial. 




Figure 23. All vertices 1 are connected through a hub in the 2nd layer. Gray thick 
lines indicat edges that have to be merged. In this manner, vertex 1 is connected to 
all other vertices. The same procedure is repeated for vertices 2 through the 3rd layer 
in order to connect vertex 2 to all the rest, and so on up to vertices n — 4. 




Figure 24. A decorated clique of 7 vertices. 

Completeness of the Ising model without magnetic fields on a 2D square lattice with 
crossings with real parameters for all Ising models with magneitc fields. In the above 
arguments we have chosen a 3D square lattice as a non-planar graph, but we can as well 
consider a 2D square lattice with crossings as our starting non-planar graph (such as 
the one in figure fT4T) . Indeed, a 3D square lattice can be constructed from a 2D square 
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lattice with crossings only by merging and deleting edges as it is shown in Example H] 
(there it is shown for the state \<p) but it also holds for the state since there also 
exist merge and deletion rules for which correspond to real parameters). Intuitively, 
they are equivalent because any edge with multiple crossings can be prepared from a 3D 
square lattice by directly merging this edge through an upper layer, as shown explicitly 
in section 15 . 11 To prepare the same edge with a 2D square lattice with crossings, 
we would merge many edges each of which has one crossing, as shown explicitly in 
Example HI Thus, we have effectively "reduced" the third dimension (or the upper 
layer) to a crossing, thereby rendering a non-planar graph which has fewer qubits and 
can prepare the same non-planar graphs. This proves that the Ising model without 
magnetic fields on a 2D square lattice with crossings is complete with real parameters 
for all Ising models. Moreover, the 2D square lattice with crossings is still polynomially 
enlarged w.r.t the target graph, growing with an overhead of (9(n 3 )(9(n 4 / 3 ) ~ 0(n 13 ^ 3 ). 

This result emphasizes the difference between planar and non-planar graphs, 
instead of 2D and 3D graphs. This difference also arises in terms of complexity, since 
the Ising model without magnetic fields on a planar graph is in the complexity class P, 
whereas the same problem on a non-planar graph is NP-hard [33]. Other non-planar 
regular graphs have been considered in the literature, such as the Kuratowskian, which 
is a 2D lattice with the crossings distributed in a sparser way than ours [33]. Notice 
that in order to transform these lattices to each other one simply needs to delete some 
crossings or merge the edges between them. 

5.2. Completeness of 3D q— state models for q— state models with pairwise interactions 

Here we show that a q— state model on a 3D square lattice is complete with real 
parameters for all q'— state models, with q' < q (both vertex and edge models) with 
pairwise interactions. The same completeness results hold for a 2q— state model on a 
2D square lattice. In order to prove the results, let us first discuss what measurements 
applied on the state \GHZ q ) (equation ( |T4l) ) correspond to real parameters. In 
particular, we shall see that the merge and deletion rules correspond to real parameters, 
whereas the creation of a crossing and the preparation of the state \xi) (cf. equation (Tl5l) ) 
require complex parameters. 

As before, we express projections in the computational basis and identify the 
coefficient of (s^) = s a ,Sb for pairwise interactions) with e /3/l ^ s<l) ' ) (see (1151) ). We 
see that only projections with all coefficients real and positive will correspond to real 
values of the Hamilton function /ij(s^). 

Merge rule. Here the merge rule consists of a certain projection on qubits belonging 
to two different GHZ states such that they become one single GHZ state. This projection 
is |00) + |11) + . . . + \q - 1, q - 1) for GHZ states of g-levels (see figure 122(a)). Note 
that it corresponds to real values of the energy /ij(s^). More precisely, it corresponds 
to assigning hi(s a ,Sb) = oo for s a = Sb, and hi(s a ,Sb) = — oo otherwise. That is, the 
coupling between the two is infinitely large when they are in the same state and is 
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infinitely small otherwise, thus forcing the particles to be in the same state and become 
one single particle. 

Deletion rule. This rule consists of deleting the interaction between two particles. 
This is achieved by the projection of two GHZ states onto (see figure 125Tb)). 

This corresponds to a vanishing coupling hi(s^) = for all spin configurations s^. The 
merge and the deletion rule are the only rules in the GHZ picture that correspond to 
real parameters. 




Figure 25. (a) Merge rule in the GHZ picture between GHZ states of q— level systems, 
(b) Deletion rule in the GHZ picture. 

Creation of a crossing. Here the creation of a crossed interaction between two pairs 
of particles is equivalent to generating two merges of interactions that cross each other 
(see figure [261 right). Since the merge and the deletion rules alone cannot generate a 
crossing, we need to start from a non-planar structure in order to obtain crossings with 
real parameters. In a system of GHZ states of q levels (see figure 126!) . a GHZ state of 
2q levels is a "non-planar" structure because the uppermost q levels can be used to 
create a crossing in the following way. For the pair of q— level GHZ states above and 
below the central one applies the merge rule going through the 2q— level GHZ state (u, 
d and c GHZ states of figure 1261) . Now the GHZ states left and right (/ and r) of the 
central one have to be merged while crossing the first one. To do this, the q— levels of 
/ are mapped to the uppermost q levels of c, i.e. the projection reads YltZo \i)M + 
Then, these upper levels of c are mapped back to the q levels of r, YH=l \l + r i)c\i)r 
(see figure 1261) . In this manner, particles I and r are effectively merged. Notice that all 
projections correspond to real Hamilton functions. That is, we have created a crossing 
corresponding to real Hamilton functions because of the extra q levels of c. Note that 
this is the minimal "non-planar" structure required to create a crossing corresponding 
to real parameters. 

Preparation of\xi)- This procedure does not correspond to a measurement rule, but 
it is a technique that we used in Secs. l3.2l and l3.3l to prove the results. In the preparation 
of q— state models, we prepared the state \xi) = Ui\0) md by constructing a circuit (in 
the one-way model) that realizes the unitary C/j. This circuit generally involves Y 
measurements as well as tilted measurements not in the range cos(6 l /2)|0) + sin(6*/2)|l), 
with 9 G [0,7r], both of which correspond to complex parameters. Therefore, we want 
to study what gates Z7j can be performed with the restriction that the coefficients of 
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Figure 26. Creation of a crossing between q— level GHZ states. It requires a non- 
planar structure if it has to be done with real parameters. The minimal "non-planar" 
structure for a q— level model is a 2q— level GHZ state. The upper and lowest GHZ 
states are merged through the lowest levels of this 2q GHZ state. The levels of the left 
particle I are mapped to the upper q levels of c, and these are mapped back to the q 
levels of r. 

all projections are real and positive (i.e. they are in the above-mentioned range). It 
seems that one cannot perform any useful gate (such as Hadamard, SWAP, CNOT,...) 
by applying projections with real and positive coefficients onto \<P2d)- m conclusion, we 
find that applying (non trivial) gates always corresponds to complex parameters. 

Now we can proceed to proving the main results of this section. 

Completeness of 3D q— state models with real parameters for all q— state models 
with pairwise interactions. In section l3~2l we mapped models with arbitrary q to an Ising 
model (q = 2) applying gates Ui that prepared \xi)- But we have just seen that these 
gates require complex parameters, and here we want to obtain completeness results with 
real parameters. Hence, we take another approach and we consider q— state models as 
candidates to be complete with real parameters for other q'— state models, with q ^ q' . 

In particular, we want to show that a q— state model on a 3D square lattice is 
complete with real parameters for all q'— state models on arbitrary graphs, with q' < q, 
with pairwise interactions. We have just seen that the merge and the deletion rule 
for q— state models in the GHZ picture correspond to real parameters. It follows that 
all observations regarding these two rules also hold for q— state models; in particular, 
one can transform a 3D square lattice made of q— level systems only with the merge 
and deletion rule into GHZ states of q— level systems with an arbitrary interaction 
pattern, but containing at most pairwise interactions. The reason why we cannot prepare 
interactions between more particles with real parameters is that their preparation seems 
to require gates which involve complex parameters. Thus, in order to prepare with real 
parameters a model with k— body interactions, the initial model must contain k— body 
interactions as well. In summary, if we allow the projections of the GHZs arranged in 
the 3D square lattice to have its maximum dimension q 2 , we can prepare with real 
parameters any other model between q— level systems involving at most two— body 
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interactions. Moreover, since these states will be prepared using merge and deletion 
rules as well, the same efficiency results hold, namely that the 3D q— state model is 
polynomially enlarged w.r.t the target model. This concludes the proof. 

Completeness of a 2D 2q— state model with real parameters for q— state models 
with pairwise interactions. Similarly as we did for Ising models, we realize that we 
require the 3D structure only because we need to create crossings by merge and deletion 
rules alone, corresponding to real parameters. However, we also know that we can use 
2q— level systems in a 2D square lattice to create create a crossing with real parameters 
between q— level systems (figure 126]) . Since we can also use the merge and deletion 
rules for this 2D 2q— level system, this proves that 2q— level system on a 2D square 
lattice is complete with real parameters for all q'— state models, with q' < q, with 
pairwise interactions. Moreover, the overhead in the system size is polynomial, because 
a polynomially enlarged 2D 2q— state model can prepare a 3D q— state model (as shown 
in Example H] for the case q — 2), and the latter needs to be polynomially enlarged w.r.t 
the target model. 

5.3. Completeness of 2D q— state vertex models for all q— state edge and vertex models 
with at most J^—body interactions 

Here we show that a q— state vertex model on a 2D square lattice (see figure [5]) with 
most general interactions in the vertices is complete with real parameters for q'— state 
vertex and edge models with at most 4— body interactions, with q' < q. Moreover, the 
2D q— state vertex model needs to be polynomially enlarged w.r.t the target model, be 
it vertex or edge model. 

In order to show that, we need to investigate what measurements can be performed 
on a vertex model in the PEPS picture that correspond to real parameters. We will 
see that the situation is different as compared to edge models: not only the merge and 
deletion rule correspond to real parameters, but also the creation of a crossing can be 
performed with real parameters. However, the performance of a gate also seems to 
require complex parameters. 

We first note that the coefficients of the projections also have to be identified with 
the coefficients of (1T51) . 

Merge rule. By applying the projection |0000) + 1 1 1 11) + . . -+\q — 1, q— 1, q— 1, q— 1) 
onto four particles belonging to different Bell pairs, we effectively merge them in the 
sense that they become a "Bell pair of four particles", that is to say, a GHZ state of 
four particles (see figure 1271(a)). Note that this measurement rule corresponds to real 
parameters. 

Deletion rule. By applying the projection onto four particles 

belonging to different Bell pairs, we effectively obtain a deletion rule, in the sense 
that these four particles become disconnected through this projection, i.e. there is 
no interaction (see figure 127Tb)). This rule also corresponds to real parameters. The 
correspondence between each of these rules and the value of the Hamilton function is the 
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same as for the merge and deletion rule for edge models in the GHZ picture (section [5.2p . 

Creation of a crossing. To create a crossing, we have to link particles up u and 
down d, and independently, right r and left I. Thus the projection ( 1 00) + |11) + . . . + 
\q — 1, q — l)) u d <8> ( 1 00) + |11) . . . + \q — 1, q — l)) r ; creates a crossing (see figure [271(c)). 
Note that this crossing is created with real parameters and this is possible because the 
projection itself is between 4 particles. 




Figure 27. Measurement rules in the PEPS picture, (a) Merge rule, (b) Deletion 
rule, (c) Creation of a crossing. The three rules correspond to real parameters. 

Performance of a gate. In section 13.31 we have proven how to map q— state vertex 
models to the 2D Ising model by performing gates Ui which prepare the states \Xi)- 
However, as we noted in section [531 the performance of gates always seems to correspond 
to complex parameters. 

Now we can proceed to prove the main results of this section. 

Completeness for vertex models. To prove that a 2D q— state model is complete 
for vertex models, we note that it can prepare other q'— state vertex models (with 
q' < q) by merging and deleting edges, and creating crossings, all of which correspond 
to real parameters. However, if the final model has an interaction between more than 4 
particles, this should be prepared using a gate, which involves complex parameters. It 
follows that in order to prepare (with real parameters) models with k— body interactions, 
the A;— body interactions must be present in the source model. This is also the reason 
why we cannot prepare models with q" > q: even though the Bell pairs with q" levels can 
be prepared, an interaction between has to be prepared via a gate, which corresponds to 
complex parameters. Hence, a 2D q— state vertex model with most general interactions 
in the vertices (i.e. of dimension g 4 ) can prepare with real parameters any q'— state 
vertex model with q' < q with at most 4— body interactions. 
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Completeness for edge models. To prove that a 2D q— state model is complete 
for edge models, we note that it can prepare a GHZ state by merging 4 Bell pairs 
(see figure 1271) . Since it can also delete and create crossings with real parameters, it 
follows that it can prepare any q'— state edge model, with q' < q, with at most 4— body 
interactions. 

We remark that in both cases the enlargement of the 2D q— state vertex model 
w.r.t the target model is polynomial, because the models are prepared using the merge, 
deletion and creation of a crossing rules, and therefore all the previous efficiency results 
hold. 

6. Relations to complexity theory 

In this work we have established reductions from general classes of classical spin systems 
to (among others) the 2D and 3D Ising model. More precisely, we have found that the 
problem of computing the partition function of a given classical spin model can be 
seen as a special instance of the problem of computing the partition function of a 2D 
Ising model with suitably tuned parameters. In many cases, we have found explicit and 
efficient reductions from these models to the 2D Ising model. 

The concept of efficiently reducing one (class of) problem(s) to another one is 
common practice in the theory of computational complexity. In fact, the computational 
complexity of certain standard spin systems, such as e.g. the Ising and Potts models, is 
well-studied. In this section we sketch (very briefly) how the results presented in this 
paper are related to some of these established complexity-theoretic results. 

The relevant complexity class in the study of partition functions is the class jfP 
("sharp-P") [34]. This class is concerned with "counting problems". More specifically, 
4fP is defined as follows: let /: {0,1}^ — > {0,1} be an arbitrary Boolean function 
(representing a decision problem) which can be evaluated in poly(iV) time; that is, 
/ represents an arbitrary problem in the class P. The complexity class #P is then 
concerned with the problem of computing the integer defined by := \{x : f(x) = 
0}|. In other words, problems in jfP are concerned with counting how many inputs 
of a given Boolean function yield the output 0. While such problems may seem rather 
innocent at first sight, they are generally very hard; in fact, there are problems in #P 
which are harder than all problems in the class NP ("non-deterministic polynomial 
time"), which itself contains numerous problems which are considered to be intractable. 
For example, the problem of counting the number of satisfying assignments to a Boolean 
3SAT formula, is a #P problem which is considered intractable (in fact, this problem 
is known to be #P-complete |34j). 

It is well known that the evaluation of the partition function of certain classes 
of classical spin systems constitutes a ^P-hard problem [35]. For example, the 
computation of the partition function of the 2D Ising model with (inhomogeneous) 
magnetic fields is known to be #P-hard [3] . This means that every problem in the class 
#P can be reduced, with polynomial computational effort, to the computation of the 
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partition function of this model. Similar results have been found for e.g. the 3D Ising 
model without magnetic fields, as well as the Potts models [35] . 

The ^P-hard completeness of the 2D Ising model hence implies that there exist 
polynomial reductions from other ^P problems (such as the evaluation of partition 
function of Ising models on arbitrary graphs, or of the Potts models) to the 2D Ising 
case. In this sense, the existence of polynomial reductions from certain spin models 
to the 2D Ising case — as we have also reported in the present work — is not a novel 
property. However, although it was indeed known that such reductions exist, making 
them explicit is by no means trivial. This is in particular illuminated by the fact that, 
for instance, the proof of #P-hardness of the 2D Ising partition function is obtained by 
reducing this problem problem to an abstract problem from graph theory (the "MAX 
CUT" counting problem), which is known to be T^P-complete. Such a proof technique 
is often used in complexity theory and acts as a kind of shortcut. Unfortunately, this 
technique does not seem to give any hint about the explicit form of the reductions from, 
say, the 3D Ising model to the 2D case. Thus, it is unclear how useful these complexity 
results are for practical purposes. 

One of the contributions of the present work is that here these reductions (e.g. to the 
2D Ising model) are made very explicit: they immediately follow from the measurement 
pattern which is needed to produce the respective stabilizer states from the (decorated) 
2D cluster state (and also this measurement pattern can be determined efficiently). It 
hence becomes possible to go beyond the abstract complexity-theoretic statement that 
the 2D Ising partition function is #P-hard, and obtain a result which may hopefully 
become useful in practice. 

Another aspect of the present results which does not seem to be covered by 
previous complexity-theoretic work, is that the reductions presented here from general 
spin systems to, for example, the 2D Ising model, are "global" in the following sense: 
given a partition function of an Ising model Za{J a b, h a ) with pairwise couplings J a b 
and local fields h a on an arbitrary graph G, we have shown that there exist a 2D 
lattice (which is polynomially larger than G) and suitable couplings J[a and h\ such 
that Z G (J ab , h a ) = Z 2 D{J'ij, hi)' (up to an irrelevant factor). In this equation, the size of 
the 2D lattice only depends on G but not on the parameters J a b and h a of the original 
model. Moreover, there is a clear relation between the original parameters {J a b, h a } and 
{JLjh'i}. In this sense, we obtain a "global" mapping from the "entire" spin model on 
G to the 2D lattice, independent of the couplings which are considered. As far as we 
are aware of, this property does not follow from, say, the #P-completeness of the 2D 
Ising model. 

7. Conclusions 

In this paper, we have utilized and extended cross-connections between classical spin 
models and measurement-based quantum computation. Using tools from quantum 
information theory, we have shown that there exist mappings between classical spin 
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systems that leave the partition function invariant. With these mappings we have been 
able to identify complete models that allow one to express the partition function of 
any (discrete) classical spin model as a special instance of the complete model. The 
mappings are in general between systems with different number of spins, that is, a 
certain enlargement of the complete model is required so that its partition function can 
specialize to other partition functions. In all relevant cases, the overhead turns out to 
be only polynomial. 

We believe that one of the main merits of our mappings and completeness results 
is that they lead to explicit constructions, which make them also interesting from the 
viewpoint of complexity theory. We have demonstrated such an explicit construction for 
the reduction of the 3D Ising model to the 2D Ising model with magnetic fields, where 
complex parameters are required. 

When allowing for complex parameters, we have found that the 2D Ising model 
with magnetic fields is complete for all classical spin models, including all q— state 
models with A;— body interactions, for any q and k (see figure [281(a) for a summary of 
the completeness results with complex parameters). As long as k and q are bounded, 
the overhead is polynomial. We have extended our analysis to vertex models, showing 
the completeness of certain 2D two— state vertex models for arbitrary models when 
allowing for complex parameters. We have also overcome the limitation of the complex 
parameters and shown that the 3D Ising model without magnetic field is complete for all 
Ising models on arbitrary graphs, using only real parameters. This implies that, for these 
mappings, a physical interpretation is possible, which means that the physical 3D Ising 
model indeed includes all other Ising models on arbitrary graphs as special instances. 
In addition, the overhead in these constructions is again only polynomial. Finally, we 
have shown the completeness of 3D q— state vertex models with real parameters for all 
q'— state models with at most 4— body interactions and with q' < q. In summary, in this 
work we have seen that, with real parameters, we can reduce models of any dimension d 
to a model with d = 3, for any fixed q and k (see figure [281(b)). However, the reduction 
to a model with d = 2, or the reduction of models with high q or k to a model with 
fixed q or k requires in our constructions the creation of crossings or the preparation of 
the state \xi), respectively, both of which correspond to complex parameters. 

Thus, one of the central questions left open in this work is whether a model can 
be found that is complete for all classical spin models, including q— state models with 
arbitrary A;— body interactions, in such a way that all parameters are real. This seems 
not to be possible for models involving only pairwise interactions, in particular we have 
not found a way to either obtain larger spin dimensions q or k— body interactions with 
k > 2 when starting from a two— state model with pairwise interactions and restricting 
ourselves to real parameters. We have, however, recently discovered such a complete 
model that involves 4— body interactions in a 3D lattice [36], where interactions take 
place between particles on the faces of cubes. 
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Figure 28. The space of all theories: the dimension d versus the number of levels 
of each particle q versus the number of many body interactions k in the theory (a) 
Completeness results with complex parameters. All theories can be reduced to the 2D 
Ising model (a model with d = q = k = 2) by invoking the universality of the cluster 
state, (b) Completeness results with real parameters. The construction of the clique 
allows to reduce models with any d to a model with d = 3, for any q and k fixed. In 
our constructions, models with higher q or k have to be reached by preparing \xi), 
which involves complex parameters. The "reductions" of models to another model 
with larger d, q or k are trivial. 
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